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Abstract — The probability density function (PDF) and cumu- 

] lative distribution function of the sum of L independent but not 
necessarily identically distributed gamma variates, applicable to 
the output statistics of maximal ratio combining (MRC) receiver 
operating over Nakagami-m fading channels or in other words to 

] the statistical analysis of the scenario where the sum of squared 
Nakagami-m distributions are user-of-interest, is presented in 
closed-form in terms of well-known Meijer's G function and 
easily computable Fox's H function for integer valued and non- 

\ integer valued diversity orders (i.e., m fading parameter of 
Nakagami fading environment), respectively. Further analysis, 
particularly on bit error rate via PDF-based approach, too is 
represented in closed form in terms of Meijer's G function and 

\ Fox's H function for integer valued diversity orders, and extended 
Fox's H function {H) for non-integer valued diversity orders. The 
proposed results complement previous elegant results that are 
either evolved in closed-form, or expressed in terms of infinite 

] sums or higher order derivatives of the diversity order m. 

I Index Terms — Gamma variates, cellular mobile radio systems, 
non-integer parameters, diversity, maximal ratio combining, 
binary modulation schemes,_bit error rate. Fox's H function, 
Meijer's G function. Fox's H function, and extended Fox's H 
function. 



I. Introduction 

IN recent times, different diversity schemes have marked an 
important impact in the arena of wireless communication 
systems. The main reason behind this is that these different 
diversity schemes allow for multiple transmission and/or re- 
ception paths for the same signal [1]. One of the optimal 
diversity combining scheme is the maximal ratio combining 
(MRC) diversity scheme where all the diversity branches are 
processed to obtain the best possible devised and improved 
single output that possibly stays above a certain specified 
threshold [l]-[3]. 

Additionally, wireless communications are driven by a com- 
plicated phenomenon known as radio-wave propagation that is 
characterized by various effects such as fading, shadowing and 
path-loss. The statistical behavior of these effects is described 
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by different models depending on the nature of the communi- 
cation environment. The wide versatility, experimental validity, 
and analytical tractability of Nakagami-m distribution [4] has 
made it a very popular fading model for performance analysis 
investigations in diversity schemes of wireless communica- 
tions (for instance, [5] among others and [6] and references 
therein). In addition, it is useful to mention that Nakagami- 
m distribution is useful to study multihop relay networks [7], 
[8]. Hence, this and such other distributions have many other 
applications in wireless communication engineering problems 
and one of those that we focus on is a communication system 
employing MRC diversity scheme undergoing this distribution 
i.e. the study of MRC diversity combining receiver operating 
over Nakagami-;7z fading channels [4], where the statistics of 
the sum of gamma random variates (RVs) or equivalently the 
sum of squared Nakagami-;7z RVs are required and moreover, 
the performance analysis of such wireless communication 
systems usually requires complicated and tedious tasks related 
to statistics as elegantly explained in details in [9]. 

The probability density function (PDF) and cumulative 
distribution function (CDF) of the sum of L independent 
but not necessarily identical (i.n.i.d.) gamma RVs has been 
elegantly investigated quite a couple of times earlier in [6], 
[9]-[ll] but the published results are sometimes given in 
rather complicated expressions in the form of single definite 
or indefinite series that renders the given expressions therein 
highly intractable or they have been derived via moment- 
generating function (MGF)-based approach. In more details, 
Moschopoulos has proposed in [11] an infinite-series repre- 
sentation for the PDF of the sum of the i.n.i.d. gamma RVs 
and Alouini et. at have extended the results of [11] in [6] 
for the case of arbitrarily correlated gamma RVs and studied 
the performance of MRC among other receivers. Commonly, 
the MGF-based approach or characteristic function (CF)-based 
approach have been followed for the performance analysis 
and derived analytical results in terms of either infinite sums 
and/or higher order derivatives of the diversity order [12]-[16]. 
This occurs as there are no simple closed-form expressions 
available in the open literature either for the PDF or the CDF 
of the sum of i.n.i.d. gamma RVs [9]. Recently, Karagiannidis 
et. al have obtained in [9] closed-form expressions for the PDF 
and the CDF of the sum of nonidentical squared Nakgami- 
m RVs or equivalently gamma RVs with integer-order fading 
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parameters but these results involve infinite series summations. 

In this work, we offer novel closed-form expressions for 
the PDF and CDF of the sum of i.n.i.d. gamma RVs or equiv- 
alently squared Nakagami-m RVs with integer-order as well 
as non-integer-order fading parameters in terms of Meijer's 
G function [17] and Fox's H function [18]-[20, App. (A.5)] 
respectively. It is noteworthy to mention that the bit error rate 
(BER) is one of the most important performance measures that 
forms the basis in designing wireless communication systems. 
Hence, we demonstrate closed-form expressions of the BER, 
as a performance metric, for binary modulation schemes, via 
PDF-based approach, of a L-branch MRC diversity receiver in 
the presence of gamma or Nakagami-m multipath fading, in 
terms of Meijer's G function and Fox's H function for integer- 
order fading parameters, and in terms of extended Fox's H 
function (H) (defined in (T.I.I) of Table I) for non-integer- 
order fading parameters. This proves the importance and the 
simplicity in the employment of those earlier derived simple 
closed-form statistical PDF and CDF expressions. These re- 
sulting easy-to-evaluate expressions also give an alternative 
form for previously known/published results obtained via CF 
or MGF-based approaches. It should be noted that all our 
newly proposed results are readily computable by Mellin- 
Barnes theorem that further corroborates the generality and 
the usefulness of the analytical frameworks introduced in this 
paper These have been checked and validated by Monte Carlo 
simulations. 

The remainder of the paper is organized as follows. Section 
II introduces the system and Section III gives novel closed- 
form expressions for the PDF and CDF of the sum of gamma 
or equivalently squared Nakagami-m RVs in terms of Meijer's 
G function and Fox's H function respectively for integer-order 
and non-integer-order fading parameters respectively. Next, 
Section IV utilizes these results presented in Section III to 
derive useful expressions for the BER, as a performance metric 
for MRC diversity receivers operating over i.n.i.d. gamma 
fading channels or equivalently Nakagami-m diversity paths in 
terms of Meijer's G function and Fox's H function for integer- 
order fading parameters, and extended Fox's H function (H) 
for non-integer-order fading parameters. And finally. Section 
V discusses the results followed by the summary of the paper 
in the last section. 

II. Nakagami-to Channel Model and Gamma 
Distribution 

A MRC based communication system with a source and 
a destination is considered with L diversity paths undergoing 
i.n.i.d. Nakagami-m fading channels as follows 

Y^aX + n, (1) 

where Y is the received signal at the receiver end, X is the 
transmitted signal, a is the channel gain, and n is the additive 
white Gaussian noise (AWGN). In a Nakagami-m multipath 
fading channel, 7 = jap follows gamma distribution. Hence, 
the channel gains experience multipath fading whose statistics 
follows gamma distribution whose PDF is given by 



^7(7) 



/TO 



exp 



7 



(2) 



where m > and > are known as fading figure 
representing the diversity order of the fading environment and 
the mean of the local power, respectively, and where r(-) 
denotes the Gamma function [21, Eq. (8.310)]. In more details, 
the parameter to quantifies the severity of multipath fading, 
in the sense that small values of to indicates severe multipath 
fading and vice versa. The instantaneous SNR of the Ith branch 
is given by 7; = [Eh /No) xi^ where xi is the envelope of the 
received signal for the lt\\ branch, Eb is the average energy 
per bit and A^o is the power spectral density of the AWGN. 

III. Closed-Form Statistical Characteristics for 
THE Sum of Gamma Random Variates 

This section presents the results on the statistical charac- 
teristics including the PDF and CDF of the sum of i.n.i.d. 
gamma variates. To best of the author's knowledge, it is 
useful to mention again that the PDF of the sum of gamma 
distributions given in the following theorem is a novel closed- 
form result not reported in the literature earlier It includes 
special cases that are used in the literature such as independent 
and identically distributed (i.i.d.) gamma RVs and/or integer- 
fading figure parameters among others [6], [9]. 

Theorem 1 (PDF of the sum of gamma or equivalently squared 
Nakagami-m RVs). Let {7;}^i be a set of i.n.i.d. gamma 
variates with parameters mi and 17;. Then, the closed-form 
PDF^ of the sum 



(3) 



for both integer-order as well as non-integer-order fading 
parameters can be expressed in terms of Fox H-function as 



1=1 





^(1) 1 


exp(y) 









(4) 



where y > 0, the coefficient sets and k G 

defined as 

k-bwcketed terms 



"(1) 



1 



TOl 



1 



TOfc 



l,TOfc , (5) 



and 



k-hracketed terms 



"(2) 



TOl 



l,TOl 



(6) 



respectively. 



'For correlated diversity branches, the statistical characteristics derivation 
and the performance analysis can be carried out in a similar fashion as 
independent fading case. For an arbitrarily correlated Nakagami-m fading 
environment, assuming that the fading parameter is common to all the 
diversity branches, the desired MGF of the sum of coiTelated gamma RVs 
can be expressed as My{s) = det(/ + sRA)^™ = Y\f=i-^iii^) = 
(1 -I- A; s)^™ where I is the L X L identity matrix. A is a positive 
definite matrix of dimension L (determined by the branch covariance matiix), 
i? is a diagonal matrix as i? = diag(ni/m, . . . , f!i/m), and A; is the 
eigenvalue of matrix RA where each eigenvalue is modeled as a gamma RV. 
Hence, on replacing 7j's with A;'s, in our presented work, we will achieve 
the results applicable to correlated diversity case. 
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TABLE I 

Representation of the Extended Fox's H function (H) and its Special Cases 



The extended Fox's H function {H) is defined as 



i,p 1 
1,9 J 



c nj=„+i{r(a, - A,s)r^ n,'=,„+i{r(i + s,^)}^ 



-z ds, 



(T.l.l) 



which contains fractional powers of F-functions. Here z may be real or complex but is not equal to zero and an empty product is interpreted as unity; C is a suitable contour, 
and positive integers p, q, m, and n satisfy the following inequalities: l<m<Q, 0<n<p, Aj > 0(j — 1, . . . , p), Bj > 0{j — 1, . . . , g) and aj (j — 1, . . . , p), 
and — 1, . . . , q) are complex parameters. The exponents aj{j — 1, . . . ,p) and bj{j — 1. . . . , q) can take on non-integer values. The poles of this integrand are 

assumed to be simple and the contour in this definition is presumed to be imaginary axis Re{s)=0 that is suitably intended in order to avoid the singularities of the gamma 
functions and to keep these singularities at appropriate sides. 



When the exponents aj 
by [lS]-[20] as 



1 for ( J — n + 1 , 



. p) and bj — 1 for {j ■ 



(A-,-Bi)i,™,(/3j,-Bj,&j)™ + i., 



.... m), the extended Fox's H function (H) reduces to the familiar Fox's H function defined 



'2^1 Jc n^_+ir(c 



(T.I. 2) 



which contains fractional powers of F-functions. Here z may be real or complex but is not equal to zero and an empty product is interpreted as unity; C is a suitable contour, 
and positive integers p. q, rn, and n satisfy the following inequalities: 1 < m < q, < n < p, Aj > 0{j — 1, ... , p), Bj > 0{j — 1, . . . , q) and (j — 1, . . . , p), 
and /3j{j — 1, .... g) are complex parameters. The exponents aj{j — 1, . . . , n) and 6j (j — m + 1, ■ ■ ■ , can take on non-integer values. The poles of this integrand 
are assumed to be simple and the contour in this definition is presumed to be imaginary axis Re(s)=0 that is suitably intended in order to avoid the singularities of the gamma 
functions and to keep these singularities at appropriate sides. 



Now, when the exponents aj 



IVj, the Fox's H function reduces to the familiar Fox //-function defined by [17] as 

n;=i r(i - a, + A,s) UT=i r(ft- - B,s) 



r I (/}i,Si),...,(/3,,B,) J nj=„ + i r(a. 



-Ai 



(T.I.3) 



which contains fractional powers of F-functions. Here z may be real or complex but is not equal to zero and an empty product is interpreted as unity; C is a suitable contour, 
and positive integers p, q, m, and n satisfy the following inequalities: l<m<Q, 0<n<p, Aj > 0{j — 1, . . . , p), Bj > 0(j — 1, . . . , g) and aj (j — 1, . . . , p), 
and /3j (j — 1, . . . , g) are complex parameters. The poles of this integrand are assumed to be simple and the contour in this definition is presumed to be imaginary axis 
Re(s)=0 that is suitably intended in order to avoid the singularities of the gamma functions and to keep these singularities at appropriate sides. 

Finally, when the exponents ctj — /3j — IVj, the Fox's H function reduces to the fainiliar Meijer's G function defined by [17] as 

1 r n"=i r(i - a, + s) n"Li r(/3, - s) 



L I /3i,...,/3, J 27riJc 



nf=„+i r(a, - s) n]=™+i r(i - ft + s) ■ 



(T.I. 4) 



which contains fractional powers of P-functions. Here z may be real or complex but is not equal to zero and an empty product is inteipreted as unity; C is a suitable contour, 
and positive integers p. q, m, and n satisfy the following inequalities: 1 < m < g, and < n < p. aj — 1, . . . , p), and — 1, . . . , g) are complex parameters. 

The poles of this integrand are assumed to be simple and the contour in this definition is presumed to be imaginary axis Re(s)=0 that is suitably intended in order to avoid the 
singularities of the gamma functions and to Iceep these singulaiities at appropiiate sides. 



Proof: In order to derive the PDF of Y, we proceed as 
follows. Firstly, the MGF 

Mj{s) =E[exp{— J s)] = / cxp(-7s)p^(7)d7 (7) 

Jo 

of a single gamma distribution is given as [1, Eq. (2.22)] 



Mj{s) = 1 + — s 



(8) 



Then, after performing some simple algebraic manipulations 
using [22, Eq. (6.1.15)], we can rewrite the MGF of a single 
gamma distribution as 



pm 1 




1 




[1 


+ n + 





(9) 



Since, 7;S are independent, the MGF of Y is the product of 
the MGF's of the 7;S, i.e. 



(10) 



Now, we express the PDF of the sum of gamma RVs, using 
the obtained MGF in (10), via inverse Laplace transform [23] 

pviy) = C-'{M{s)} = </ MY{s)eMs)ds (11) 
27ri Jc 

that produces a Mellin-Barnes contour integral [24] represen- 
tation as (similar to (T.l.l), Table I) 



1=1 
1 

27ri 



cxp(s)(is. (12) 



1=1 



Hence we use this obtained result and perform some simple 
rearrangements on the r(.) terms in the Mellin-Barnes contour 
integral representation to express the closed-form PDF of the 
sum of gamma RVs Y, valid for both integer-order as well 
as non-integer-order fading parameters, in terms of Fox's H 
function as given in (4). ■ 
The motivation and possibility that led to the above result 
presented in Theorem 1 was the representation of the PDF a 
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single gamma RV in terms of Meijer's G function as discussed 
below in Corollary 1. 

Corollary 1 (PDF of a single gamma RV). Let {7;} be any 
i.n.i.d. gamma variate with parameters mi and fJ;. Then, the 
closed-form PDF of this single gamma RV can be expressed 
as 



Pii ill) 



mi 



G 



m;,0 



exp(-7;) 



(1) 

mi 
(2) 



(13) 



where the coefficient sets and ^^^\ G N are defined as 



k'times 



and 



k-times 



(2) _ / ruk 



(14) 



(15) 



respectively. 



Proof: The PDF of a single gamma RV can be expressed 
using the obtained MGF in (9), via inverse Laplace transform 
[23] that produces a Mellin-Barnes integral [24] representa- 
tion. Hence we use this obtained result to express the PDF of 
a single gamma RV in an alternative form, in terms of Meijer's 
G function, as expressed in (13) ■ 

To our best knowledge, the Fox's H function is not available 
in any standard mathematical packages. As such, we offer 
in Table II an efficient Mathematica® implementation of 
this function (similar to [7], [25], [26]) in order to give 
numerical results based on (4). With this implementation, 
the Fox's H function can be evaluated fast and accurately. 
This computability, therefore, has been utilized for different 
scenarios and is employed to discuss the results in comparison 
to respective Monte Carlo simulation outcomes. 

Fig. 1 and Fig. 2 present the PDF and the logarithmic PDF 
respectively of the output SNR obtained from the exact closed- 
form expression (4) and show a perfect match between this 
obtained closed-form analytical result and the one obtained via 
Monte Carlo simulations for varying L's (i.e. L — 3,4, 5), and 
their respective fixed fading parameters mi — 0.6, TO2 ~ 1.1, 
ma — 2, m4 — 3.4, and m^ = 4.5. 

For additional verification purposes, we simplified the PDF 
expression presented in [9, Theorem 1, Eq. (6)] for the sum of 
squared Nakagami-;7z or equivalently gamma RVs, for a special 
case with L — 2 i.e. a dual-branch MRC diversity combining 
receiver based wireless communication system undergoing 
Rayleigh fading i.e. mi = 1 and 7712 = 1. We obtained the 
following simplified PDF expression. 



1 








cxp ( - — 











(16) 



where, rj^s are equivalent to our Q[s and z is the RV as 
opposed to our RV y. This obtained expression, when plotted 
against our results, with similar specific values, found a perfect 
match along with the Monte Carlo simulations and hence 
further conforming our results as shown in Fig. 3. 



TABLE II 

Mathematica® Implementation of Fox's H function 



(*Fox H- Bar -Function Implementation*) 

Clear [x, Q] ; 

(♦Exception*) 

FoxHBar : : InconsistentCoef fs = "Inconsistent i 



FoxHBar [a_, b_, z_ 



= Module 



{Z, 3, Pa, Pb, Qa, Qb, M, R, value}, 
{*6aimna product terms*) 
Pa = Function [u. Product [ 

Power [Gamma [1 - a[ [1, n, 1]] + ua[[l, n, 2]]], a[[l, n, 3]]], {n, 1, Length [a [ [1] ; 
Qa = Function [u. Product [Gamma [a [ [2 , n, 1]] - ua[[2, n, 2]]], {n, 1, Length [a [ [2] ] 
Pb = Function [u. Product [Gamma [b [ [1, n, 1]] - ub[[l, n, 2]]], {n, 1, Length [b[ [1] ] ; 
Qb = Function [u. Product [ 

Power[Gamma[l-b[[2, n, 1]] +ub[[2, n, 2]]], b[[2, n, 3]]], {n, 1, Length [b[[2] 
M = Function [u, Pa[u] Pb[u] /Qa[u] /Qb[u]] ; 
{•Contour Limiter*) 

(*Depends on numerator argument i.e. it 

must be at least half of the least valued gamma arguments*) 
R = -1; 

(♦Assignment and Declaration*) 
Z = z; 

(*Final Evaluation*) 
1 

value = 

2 TTl 



NIntegrate [M[s] Z% {s, -50 - I 100, 
(♦Returning back the value*) 
Return [value] ; 

{*End of FoxHBar*} 



R-IlOO, R + IlOO, -50 + 1 100), MaxRecursion -» 55] , 



Comparison between Analytical and Simulation Results 



Simulation 
Analytical 




Output Signal-to-Noise Ratio y 



Fig. 1 . Comparison between PDFs obtained analytically and via Monte Carlo 
simulations for varying branches L and respective fixed fading parameters for 
these channels as mi = 0.6, m2 = 1.1, = 2, = 3.4, and = 4.5. 



Corollary 2 (PDF of the sum of gamma or equivalently 
squared Nakagami-m RVs for integer-order fading parame- 
ters). It is worth mentioning that the closed-form expression in 
(4) simplifies to the following expression (18) /or integer-order 
fading parameters via simple algebraic manipulations. 



py{y) = n 



^^1 



mi 



cxp(-?;) 



(1) 

K 

(2) 



(18) 
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fc-bracketed terms 



m 1 -times m -times 



Comparison between Analylical and Simulation Resuits 




5 10 15 

Output Signai-lo-Noise Ratio y 



Fig. 2. Comparison between PDFs obtained analytically and via Monte Carlo 
simulations, on log scale, for varying branches L and respective fixed fading 
parameters for these channels as mi = 0.6, m2 = 1.1, ma = 2, m4 = 3.4, 
and m5 = 4.5. 



Comparison between Analytical and Simulation Resuits 

0.4 I , , , , , , 




01 23456789 10 

Output SIgnal-lo-Nolse Ratio y 



Fig. 3. Comparison between PDFs obtained analytically, via Monte Carlo 
simulations, and via [9, Theorem 1, Eq. (6)] with L = 2 and fading parameters 
for these channels as mi = 1 and m2 = 1 i.e. Rayleigh fading channels. 



where. 



L 

1=1 

(1) 



(19) 



is an integer, the coefficient set , k ^ N is as defined in 

(2) 

(17) and the coefficient set , k is defined as 



k-bmcketed terms 



(2) 



mi 



mi 



( nm. 



( VIE. 



(20) 

where K is the total number of gamma or equivalently squared 
Nakagami-m RVs i.e. L number of total branches for our 
specific wireless communication system being considered here. 

Now let us consider some special cases in order to check 
the correctness and accuracy of (18). These special cases give 
a further insight to the above obtained results and assist in 
understanding the rest of the results presented in this work. 

Special Case 1 (Sum of two exponential RVs). Let us assume 
that we have two i.n.i.d. gamma RVs with fading figures nii ~ 
1 and 771,2 = 1 and average powers fti and f72. Substituting 



these parameters in (18) results in 



Pviy) 



1 



-G 



2,0 



cxp(-2/) 



1 



1 

1 1 



1 

02 



(21) 



Then, using the Meijer's G identity given in [20, Eq. (1.142)] 
and then using [27, Eq. (07.23.03.0227.01)], (21) readily 
reduces to [28, Sec. 5.2.4] 



Pviy) 



exp(. 



-JL) 



(22) 



Special Case 2 (Sum of L exponential RVs). Let us assume 
that we have L i.n.i.d. gamma RVs with fading figures mi = 1 
for all ? G {1, 2, 3, . . . , L} and average powers ili ^ ilk for all 
k,l € {1, 2, 3, . . . , i}. Substituting these parameters in (18) 
results in 



Pviy) = 



1 



nti 



-G 



L,0 
L,L 



exp{-y) 



1 



1 

Sli ' ■ 

1 



,1 + 



(23) 

Then, performing some algebraic manipulations using the 
Meijer's G identity given in [27, Eq. (07.34.26.0004.01)] and 
[27, 07.31.06.0017.01], we simplify (23) to [28, Eq. (5.8)] 



Pviy) 




(24) 
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It is appropriate to mention that the steps followed in 
Special Case 1 and Special Case 2 can also be readily used 
to obtain the CDF for the same special cases. 

Theorem 2 (CDF of the sum of gamma or equivalently 
squared Nakagami-m RVs). The CDF of Y for both integer- 
order as well as non-integer-order fading parameters can be 
closely expressed in terms of Fox H-function as 

L 



Comparison between Analytical and Simulalion Results 



n 



mi 



H 



L+1,L+1 



exp(?/) 



(0,1,1) 



where the coefficient sets S^, 
(5) and (6) respectively. 



"(2) 

and SI are defined earlier in 



(25) 



Proof: In order to derive the CDF of Y , we proceed as 
follows. 

We integrate the PDF expressed in (13) from through 
7 and obtain the CDF for a single gamma RV, in terms of 
Meijer's G function, as 



^7(7) 



G 



m+1,0 
m+l,m+l 



exp(-7) 



$W,1 



(26) 



where the coefficient sets and are defined earlier in 
(14) and (15) respectively. 

Now, performing a similar integral operation on (4), uti- 
lizing a similar explanation as presented in the proof of the 
PDF of the sum of gamma RVs i.e. Theorem 1 to obtain (4) 
from (10), and further making some simple modifications to 
the Mellin-Barnes integral representation to satisfy the exact 
definition of the Fox's H function, we obtain a final closed- 
form result for the CDF of Y , valid for both integer-order as 
well as non-integer-order fading parameters, in terms of Fox's 
H function as presented in (25). ■ 

Therefore, Fig. 4 and Fig. 5 present the CDF and the 
logarithmic CDF respectively of the output SNR obtained 
from the exact closed-form expression (25) and show a perfect 
match between this obtained closed-form analytical result and 
the one obtained via Monte Carlo simulations for varying L' s 
(i.e. L = 3,4, 5), and their respective fixed fading parameters 
mi = 0.6, 1712 = 1.1, TO3 = 2, 1714 = 3.4, and 771,5 = 4.5. The 
logarithmic plots were selected to display the accuracy of the 
matched results. 

Corollary 3 (CDF of the sum of gamma or equivalently 
squared Nakagami-m RVs for integer-order fading parame- 
ters). It is worth mentioning that via simple algebraic ma- 
nipulations, the expression in (25) simplifies to 



1=1 



pl + K,0 



exp(-2/) 



(27) 



for integer-order fading parameters where the coefficient sets 



'^^fP and '^)^' are defined earlier in (17) and (20) respectively. 

Moreover, when the probability that the instantaneous MRC 
output SNR falls below a given threshold ?/th, we encounter 
a situation labeled as outage and it is an important feature 
to study outage probability (OP) of a system. Hence, another 
important fact worth stating here is that the expression derived 



,(2) 




Output Signal-to-Noise Ratio y 



Fig. 4. Comparison between CDFs obtained analytically and via Monte Carlo 
simulations for varying branches L and respective fixed fading parameters for 
these channels as rrii = 0.6, m2 = 1.1, "13 = 2, m4 = 3.4, and = 4.5. 



Comparison between Anaiytical and Simulation Resuits 




Simulation 
- Analytical 



Output Signal-to-Noise Ratio y 



Fig. 5. Comparison between CDFs obtained analytically and via Monte Carlo 
simulations, on log scale, for varying branches L and respective fixed fading 
parameters for these channels as mi = 0.6, m2 = 1.1, ms = 2, mi = 3.4, 
and m5 = 4.5. 



in the Corollary 3 also serves the purpose for the expression 
of OP of MRC diversity combining receivers based wireless 
communication system that is experiencing i.n.i.d. Nakagami- 
m fading channels or in other words, when the desired user is 
subject to Nakagami-m fading, the probability that the SNR 
falls below a predetermined protection ratio yth can be simply 
expressed, for both integer-order as well as non-integer-order 
fading parameters, by replacing y with yt\^ in (25) as 

Pout(2;th) = Py(yth)- (28) 
Employing similar substitutions, all the other respective ex- 
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TABLE III 

Conditional Error Probability (CEP) Parameters 



Modulation 


V 


1 


Binary Frequency Shift Keying (BFSK) 


0.5 


0.5 


Binary Pliase Shift Keying (BPSK) 


0.5 


1 


Differential Phase Shift Keying (DPSK) 


1 


1 



pressions of CDF can be utilized for OP such as replacing 7 
with 7th in (26) and/or replacing y with yth in (27). 

IV. Applications to the Performance of Diversity 
Combining Receiver Systems 

This section presents the results on the performance analy- 
sis, in particular BER analysis, of the sum of i.n.i.d. gamma 
variates or equivalently squared Nakagami-m variates. In here, 
we limit ourselves to binary modulation schemes. 

In MRC combining scheme, all the branches are selected at 
the output. In our case, for a L-branch MRC diversity receiver, 
the signal-to-noise ratio (SNR) y, is given by 



2/ = 71 H I-7L. 



(29) 



The most straightforward approach to obtain BER Pe for 
MRC is to average the conditional error probability (CEP) 
Pe (e|y) for the given SNR given by 

r(p,gy) 



over the PDF of the combiner output SNR [1] i.e. 



P. = 



Pe {Av)Py [y] dy. 



(30) 



(31) 



The expression in (30) is a unified CEP expression for coherent 
and non-coherent binary modulation schemes over an AWGN 
channel [29]. r(-, ■) is the complementary incomplete gamma 
function [21, Eq. (8.350.2)]. The parameters p and q in (30) 
account for different modulation schemes. For an extensive list 
of modulation schemes represented by these parameters, one 
may look into [30] or refer to Table III. 

Theorem 3 (BER of a L-branch MRC system operating over 
Nakagami-m fading channels for binary modulation schemes). 
The BER of a L-branch MRC diversity combining receiver 
wireless communication system running over Nakagami-m 
fading channels, valid for both integer-order as well as non- 
integer-order fading parameters and for any binary modu- 
lation scheme including BFSK, BPSK, and DPSK, can be 
expressed in closed-form, in terms of extended Fox H-function 
(H), as 



P. = 



where 



and 



1=1 



H 



{l-q,l,p)Xl 

C2, {-q, 



Ti'\ (1,1,1), 
(0,l,l),Ti'\ 



(32) 



(33) 



and where the coefficient sets T^j^^ and t[,^^, k £ N are defined 



as 



k'bracketed terms 



T 



and 



(1) 



mi 



1 + — ,l,mi + — ,l,mfc , (35) 



k'bracketed terms 



(36) 



respectively. 



Proof: Utilizing (31) by substituting (30) and (4) into it 
and performing some simple manipulations along with some 
simple rearrangements of r(.) function terms, we get an exact 
closed-form result of the integral valid for both integer-order as 
well as non-integer-order fading parameters and for any binary 
modulation scheme including BFSK, BPSK, and DPSK, in 
terms of extended Fox's H function {H), as presented above 
in (32), Theorem 3. ■ 

Corollary 4 (BER of a L-branch MRC system operating 
with Nakagami-;7z integer-order fading channels for BFSK 
and BPSK binary modulation schemes). The above presented 
BER expression in (32), Theorem 3 is simplified, via simple 
algebraic manipulations, to the following closed-form expres- 
sion when considering BFSK and BPSK binary modulation 
schemes with only integer-order fading parameters. It is rep- 
resented in terms of Fox H-function as 



P. = 



where 



and 



9 11 I f], 



1=1 



H 



K+1,1 

K+2,K + 2 



(1 - 9, 1,P),X1 

X2,(-g, i,p) 



Xi=AW,(1,1), 
X2 = (0,1),A(2), 



(37) 



(38) 



(39) 



and where the coefficient sets A^^' and A{.^', fc G N are 
defined in (40) and (41) respectively. K is the total number of 
gamma or equivalently squared Nakagami-m RVs i.e. L num- 
ber of total branches for our specific wireless communication 
system being considered here. 

Corollary 5 (BER of a L-branch MRC system operating over 
Nakagami-m integer-order fading channels for DPSK binary 
modulation schemes). Further down the line, the BER expres- 
sion in (37) is simplified, via simple algebraic manipulations, 
to the following closed-form expression when considering only 
DPSK binary modulation scheme for integer-order only fading 
parameters. This expression is represented in terms of Meijer 
G-function as 



1=1 



mi 



G 



K+P+l.ft+P+l 



(2) ^(2) 



(42) 

where the coefficient sets ^''^^ and ^''jf^ are defined earlier in 
(34) (17) and (20) respectively, the coefficient sets o[3'' and o^^\ 
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fc-bracketed terms 



m 1 -times m -times 



fc-bracketed terms 



m 1 -times m k -times 




k gN are defined as 

k'times 

0^ ='(l-g),.?.,(l-g); (43) 

k~times 

of = (44) 

respectively. 

V. Results and Discussion 

The numerical results for BER of MRC diversity combining 
receiver scheme with L-diversity over i.n.i.d. gamma or equiv- 
alently squared Nakagami-m fading channels are presented in 
this section. 

To our best knowledge, the extended Fox's H function 
{H) is not available in any standard mathematical packages. 
As such, we offer in Table Table IV an efficient Mathe- 
matica® implementation of this function in order to give 
numerical results based on (32), (37), and/or (42). With this 
implementation, the extended Fox's H function {H) can be 
evaluated fast and accurately. This computability, therefore, 
has been utilized for different digital modulation schemes and 
is employed to discuss the results in comparison to respective 
Monte Carlo simulation outcomes. 

The average SNR per bit in all the scenarios discussed is 
assumed to be equal. In addition, different digital modulation 
schemes are represented based on the values of p and q 
where p = 0.5 and q = 1 represents BPSK, p = 1 and 
q = 1 represents DPSK, and BFSK is represented by p = 0.5 
and q = 0.5. In Monte Carlo simulations, the gamma or 
equivalently squared Nakagami-m fading channel generation 
is readily available in MATLAB. 

We observe from Fig. 6 that this implemented computability 
of extended Fox's H function (H) provides a perfect match to 
the MATLAB simulated results and the results are as expected 
i.e. the BER decreases as the signal-to-noise ratio (SNR) 
increases. Its important to note here that these values for the 
parameters were selected randomly to prove the validity of 
the obtained results and hence specific values based on the 
standards can be used to obtain the required results. 

Furthermore, it can be seen from Fig. 6 that, as expected, 
BPSK outperforms the other modulation schemes and DPSK 
outperforms BFSK. Similar results for any other values of m's 
can be observed for the exact closed-form BER for L-diversity 



TABLE IV 

Mathematica® Implementation of the Extended Fox's H 
function (h) 



(^Extended Fox H-Bar-Function*) 
Clear [x, Q] ; 

ExtendedFoxHBar : : InconsistentCoeff s = " Inconsistent coefficients 1 " ; 
ExtendedFoxHBar [a_, b_, a ] :- Module |^ 

{Z, s. Pa, Pb, Qa, Qb, M, R, Rmax, Rmin, value}, 
(*Gaiiinia product terms*) 
Pa - Function [u. Product [ 

Power[Gamina[l-a[[l, n, 1] ] - u a [ [1, n, 2] ] ] , a [ [1 , n, 3 ] ] ] , {n, 1 , Length [a [ [1] ]])]] ; 
Qa ^ Function [u. Product [Power [Gamma [a [ [2,n, 1]] + ua[[2,n, 2]]],a[[2,n, 3]]], 

{n, 1, Length[a[[2]]]}]]; 
Pb ^ Function[u, Product [Power [Gamma [b [ [1, n, 1]] + ub[[l, n, 2]]], b[[l, n, 3]]], 

{n, 1, Length[b[[l]]]}]]; 
Qb ^ Function [u. Product [Power [Gamma [1 - b[[2, n, 1]] - ub[[2, n, 2]]], b[[2, n, 3]]], 

{a, 1, Length [b[[2]]]}]]; 
M = Function [u. Pa [u] Pb [u] / Qa [u] / Qb [u] ] ; 
(itContour Limiter*) 

(^Depends on numerator argument i.e. it 

must be at least half of the least valued gamma arguments*) 
{Rmin, Rmax) ^ {Max [-Min [b [ [1, All, 1]] /b[[l, All, 2]]], -Infinity], 

Min[Min[ (1 - a[ [1, All, 1] ] ) / a[[l. All, 2] ] ] , Infinity] }; 
If [Rmin ~ - Infinity Rmax ^ Infinity, Rmin - Rmax - 1] ; 
If [Rmin ^ - Infinity Rmax =: Infinity, Rmax - Rmin + 0.1] ; 
If [Rmin -- Rmax, Message [ExteadedFoxHBar : : InconsistentCoef f s ] ] ; 
R = Mean [ { Rmax , Rmin } ] ; 
(^Assignment and Declaration*) 
Z ^ z; 

(*Final Evaluation*) 

1 

KIntegrate [M[s] Z", {s, -50 - I 50, R-I50, R + I 50, -50 + 150), MaxRecursion ^ 55] ; 

2 TTI 

(^Returning back the value*) 
Return [value] ; 

]■■ 

(*End of ExtendedFoxHBar*) 



i.n.i.d. gamma or equivalently squared Nakagami-m channels 
presented in this work. 

VI. Concluding Remarks 

We derived novel closed-form expressions for the PDF and 
the CDF of the sum of i.n.i.d. gamma or equivalently squared 
Nakagami-m RVs in the case of both integer-order as well 
as non-integer-order fading figure parameters. An interesting 
finding is that these expressions can be written in terms of 
special functions, specifically Fox's H function distributions. 
Based on these statistical formulas obtained, and following the 
exact PDF-based approach, we analyzed the performance of a 
MRC diversity combining receiver based wireless communica- 
tion system operating over i.n.i.d. Nakagami-fading channels 
and important performance metrics such as OP and BER were 
expressed in closed form and hence this serves as the key 
feature along side the novel statistical derivations of PDF and 
CDF 

An exact closed-form expression for the BER performance 
of different binary modulations with L-branch MRC scheme 
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Comparison between Analyiical and Simulation Results 
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Fig. 6. Average BER of BPSK over i.n.i.d. gamma or equivalently squared 
Nakagami-»! fading channels with L = 5-branch MRC and fading pai'ameters 
for these channels as m,\ = 0.6, r?i2 = 1.1, ma = 2, 7714 = 3.4, and 
ms = 4.5. 



over i.n.i.d. gamma or equivalently squared Nakagami-m fad- 
ing channels was derived. The analytical calculations were 
done utilizing a general class of special functions including 
Meijer's G function. Fox's H function, and extended Fox's 
H function {H). In addition, this work presents numerical 
examples to illustrate the mathematical formulation developed 
in this work and to show the effect of the fading severity and 
unbalance on the system performance. 

Our results improve on previously published results that are 
either in the form of infinite sums or higher order derivatives 
of the fading parameter or recursive solutions. This approach 
provides a novel closed-form expression as compared to earlier 
work in the literature. Having the computability for these 
functions such as Fox //-function and extended Fox's H 
function (//) and/or the decomposition of these functions into 
univariate Meijer's G function, this approach will form the 
basis for further analysis in similar proceedings. 
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